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are calculated within the lattice-constrained dispersion quark model: the 
form factors are given by dispersion representations through the wave 
functions of the initial and final mesons, and these wave functions are 
chosen such that the B K* transition form factors agree with the 
lattice results at large q^. We calculate branching ratios of semilep- 
tonic B K, K* transition modes and study the sensitivity of ob- 
servables to the long-distance contributions. The shape of the forward- 
backward asymmetry and the longitudinal lepton polarization asymmetry 
are found to be independent of the long-distance effects and mainly de- 
termined by the values of the Wilson coefficients in the Standard Model. 
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I. INTRODUCTION 

The investigation of rare semileptonic decays of the B meson induced by the flavour-changing neutral current 
transitions h s represents an important test of the Standard Model (SM) and its possible extensions. Rare 
decays are forbidden at tree level and occur at the lowest order only through one- loop diagrams. This fact 
opens the possibility to probe at comparatively low energies the structure of the electroweak theory at large 
mass scales, thanks to the contributions of virtual particles in the loops. Moreover, rare b ^ s transitions are 
expected to be sensitive to possibile new interactions, like those provided, e.g., by supersymmetric theories, two 
Higgs-doublet, top-color and left-right models. These interactions govern the structure of the operators and the 
corresponding Wilson coefRcients, which appear in the AB — 1 effective electroweak Hamiltonian describing 
the b —^ s transitions at low energies. 

A recent experimental observation of exclusive |Q] and inclusive ||] radiative decays, B K*j and B 
Xs"f, have prompted a lot of theoretical investigation on rare semileptonic B decays. However, in the case 
of exclusive decays any reliable extraction of the perturbative (short-distance) effects encoded in the Wilson 
coefficients of the effective Hamiltonian ||-|^ requires an accurate separation of the nonperturbative (long- 
distance) contributions, which therefore should be known with high accuracy. The theoretical investigation 
of these contributions encounters the problem of describing the hadron structure, which provides the main 
uncertainty in the predictions of exclusive rare decays. 

In exclusive B ^ K, K* decays the long-distance effects in the meson transition amplitude of the effective 
Hamiltonian are encoded in the meson transition form factors of bilinear quark currents. Various theoretical 
frameworks have been applied to the description of meson transition form factors; among them we should 
mention constituent quark models QCD sum rules |l2| - [l4[ , lattice QCD ||l5|-|l7[|, approaches based on 

the heavy-quark symmetry [ p^ and analytical constraints |ig| ]. 

Lattice QCD simulations, because of its most direct connection with QCD, are expected to provide the most 
reliable results. Although it is not possible to place the 6 quark directly on the lattice, a constrained extrapolation 
in the heavy quark mass jlTt allows to determine reliably the form factors for B decays. A present limitation is 
that lattice calulations do not yet provide the form factors in the whole accessible kinematical decay region: the 
daughter light quark produced in b decay cannot move fast enough on the lattice and one is therefore limited 
to the region of not very large recoils. For obtaining form factors in the whole kinematical decay region one 
can use extrapolation procedures based on some parametrizations of the form factors. For instance, in [ pTf a 
simple lattice-constrained parametrization based on the constituent quark picture and pole dominance is 
proposed. Anyway, a reliable knowledge of form factors in some region is already a substantial step forward, 
which provides firm constraints for the results of other approaches. 

QCD sum rules give complementary information on the form factors as they can calculate the latter at not 
very large momentum transfers. However in practice various versions of the QCD sum rules give remarkably 
different predictions, being strongly dependent on the technical subtleties of the particular version. A recent 
analysis disregards the three-point sum rules in favour of the light-cone sum rules. On the other hand, 
the light-cone sum rules involve more phenomenological inputs and the results turn out to be sensitive to the 
particular distribution amplitude of the light meson used in the evaluation of the sum rule and to the model 
adopted for the subtraction of the continuum (cf. [|l3| and @). 
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Constituent quark models (QM) have proved to be a fruitful plienomenological method for the description of 
heavy meson transitions. An attractive feature of the approaches based on the concept of constituent quarks is 
the suggestion of a simple physical picture of the decay process, based on the following phenomena responsible 
for the soft physics: i) the chiral symmetry breaking in the low-energy region generating the constituent quarks; 
ii) a strong peaking of the soft (nonperturbative) hadronic wave functions in terms of the relative constituent 
momenta with a width of the order of the confinement scale; and iii) the dominance of the contribution of the 
Fock state components with the minimal number of constituents, i.e. qq component in mesons. An important 
shortcoming of the quark model predictions for the form factors is a strong dependence of the results on the 
QM parameters 

Thus we can see that none of the above-mentioned approaches can provide accurate form factors in the whole 
kinematically accessible region of B decays. In this situation a combination of the results of different approaches 
can be efficient for obtaining reliable predictions. 

Our approach to the calculation of the B —> K, K* transitions is based on the dispersion quark model 
pl] , |2^ . The transition form factors are given by relativistic double spectral representations through the wave 
functions of the initial and final mesons both in the scattering and the decay regions. The form factors of the 
dispersion quark model develop the correct heavy-quark expansion at leading and next-to-leading 1 /toq orders 
in accordance with QCD for the transitions between heavy quarks p3| , p4| . For the heavy-to- light transition the 
form factors of the dispersion quark model satisfy the relations between the form factors of vector, axial-vector, 
and tensor currents valid at small recoil |^. Thus the form factors of the dispersion quark model obey all 
known rigorous theoretical constraints. A possibility to calculate directly the form factors in all the decay region, 
avoiding in this way any extrapolation, is an important advantage of this formulation of the quark model. The 
main results of this work are as follows: 

• We present a dispersion quark model calculation of the B K, K* transition form factors in the whole 
kinematical range of q^ . Adopting the quark masses and the wave functions of the Godfrey- Isgur model 

for the hadron spectrum with a switched-off one-gluon exchange (OGE) potential for taking into 
account only the impact of the confinement scale, we have found that the resulting form factors are in 
good agreement with the lattice simulations at large q^ . Thus we expect to provide reliable form factors 
in the whole decay region. 

The dispersion quark model form factors for the B — > K* transition agree favorably in the whole range of 
< 9^ < {Mb — M^)"^ with a lattice-constrained fit based on the constituent quark picture ||] and an 
assumption on a single-pole behavior of the form factor Ai{q'^). On the other hand the parametrizations 
based on heavy-quark symmetry (HQS) also give reasonable results if one assumes the leading-order 
expressions for the form factors and replaces the universal process-independent Isgur- Wise (IW) function 
with process-dependent form factors, Cb^k and Cb^k*, related to the B —f K and B — > K* transitions, 
respectively. The latter are found to differ strongly from each other and from the asymptotic IW function. 

An important consequence is that both our QM calculations of the form factors and the lattice-constrained 
parametrization of Ref. |p7[| as well as the use of the heavy-quark symmetry relations between the form 
factors predict a quite similar behavior for the forward-backward and lepton polarization asymmetries. 

• We derive formulas for the differential decay rates and asymmetries in exclusive rare semileptonic decays 
of heavy mesons for the case of massive leptons and taking into account a nonzero mass of the daughter 
quark produced in the rare b transition. For massless leptons and/or in the limit rn^ — > our formulas 
reproduce known results. 

• We present a detailed analysis of non-resonant decay rates and asymmetries in — > {K,K*) [t'^i' ,vv) 
decays within the SM adopting our QM transition form factors. For comparison we also perform calcu- 
lations with the lattice-constrained form factors of Ref. . The decay rates evaluated in both models 
are found to be in agreement with each other, while the differential dilepton distributions are sensitive to 
subtle details of the g^-dependence of the transition form factors. It is found that the lepton polarization 
asymmetry (Pl) as well as the shape of the forward-backward asymmetry {Apb) are largely independent 
of the long-distance contributions and determined only by the values of the Wilson coefficients in the 
SM. Such features make both ApB and Pl good candidates for testing the Standard Model and probing 
possible New Physics. 

The paper is organized as follows. In Section II the SM operator basis, describing the b — > s£^£~ and b svv 
transitions, is briefly presented. In Section III the meson transition form factors are considered. Section IV 
presents the differential rates, lepton spectra and lepton asymmetries for the rare B {K,K*) {£~^£~,i'D) 
decays including the case of massive leptons. Section V gives numerical analysis of the lepton spectra and 
lepton asymmetries in exclusive rare i3-meson decays in the SM. Conclusion summarizes the results and gives 
an outlook. 
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II. THE OPERATOR BASIS 



The effective weak Hamiltonian, which describes the b — > s£^£ transition, has the foUowing form ||] 

Gp 



(1) 



where Gp is the universal Fermi constant, the quantities Ci(/i) are the Wilson coefficients, obtained after 
integrating out the heavy particles, and the O^'s are the basis operators; the sign of the Wilson coefficients 
Gi{iJ,) is determined as in the work C2{Mw) = —1- Within the SM, the operators providing the main 
contribution to rare decays are 

01 = (Sq7''(1 - 75)&a) (C/37^(1 - 75)c/3) , 

02 = (Sa7^(l - 75)^/3) (c/37p(l ~ 75)Ca) , 

g 

Or-y = ■^-^Sa<J^l,[mb{^i){l +75) + ms(^)(l - 75)]&q i^'"', 

OTT 

^2 



„2 



r(Sa7^(l-75)&a)^7A.^ 



(2) 



In Eq. the renormalization scale fi is usually chosen to be ^ ~ mb in order to avoid large logarithms in the 
matrix elements of the operators Oi . The Wilson coefficients Ci reflect the speciflc features of the theory at large 
mass scales; they are calculated at the scale fi ~ Mw and then evolved down to /i ~ by the renormalization 
group equations. The analytic expressions for Ci{fi) in the SM can be found, e.g., in Q|. In what follows, the 
values of the Wilson coefficients at the scale fi ~ nib ~ 5 GeV are |^,^: Ci(mf,) = 0.241, C2(mb) = —1.1, 
C7j{mb) = 0.312, Cgvimt) = -4.21 and Cioa(to6) = 4.64. 

The four-quark operators Oi and O2 generate both short- and long-distance contributions to the effective 
weak Hamiltonian (ij). Both contributions can be taken into account by replacing CgvinT-b) with an effective 
coefficient Gg{/ {mb,q'^) given by [D 



Cl(,^{mb, q^) = Gw{mb) + [3Ci(mb) + C2(mfc)] • 



TOc q_ 

rub ' ml 



E 



-iMv-T 



V^^ Vi 



(3) 



where is the invariant mass squared of the lepton pair. The short-distance contributions are contained in 
the function h(rnc/mb,q^ /rn^), which describes the one-loop matrix element of the four-quark operators Oi 
and O2 (see, e.g., Q for its explicit expression). The long-distance contribution, related to the formation of 
intermediate cc bound states, is usually estimated by combining the factorization hypothesis and the Vector 
Meson Dominance assumption P||] ; phenomenological analyses B suggest that in order to reproduce correctly 
the branching ratio BR(B ^ J/?/'X e+e'X) = BR{B J/i'X) ■ BR(J/V' ^ the fudge factor k, 

which appears in Eq. (^) to correct phenomenologically for inadequacies of the factorization -I- VMD framework, 
should satisfy the approximate relation: k [3Ci(TOh) + C2(mf,)] w 1. To sum up, the effective weak Hamiltonian 
has the following structure (cf. UJ|,@1) 



Heffib^ sl+r) 



Gp OLf-rn 
72^ 



((mfc + ms)isia^^q''b) + {mb - ms){sia^^q''j5b)) (Ij'^l) (4) 



+ Cli^'imb, q'){sj^,{l - J^miri) + aioA(mfc) (57^(1 - 75)fo)(^7^750 



Heffib^ sw) 



Gp Q!e 



-VtbV:,X{^t){-s^-1^{l - 75)fc")(i^7'^(l - 75)^^) 



72 27r sin^ dw 

where Xt — {mt/MwY and X{xt) is given in 0. At — 176 GeV one has X{xt) = 2.02. 



(5) 
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III. MESON TRANSITION FORM FACTORS. 



The long-distance contribution to _B {K,K*) decays is contained in the meson matrix elements of the 
bilinear quark currents of Ti-eff, i-e. in the relativistic invariant transition form factors of the vector, axial- 
vector and tensor currents. In rare semileptonic decays there is another long-distance effect, known as the 
weak annihilation, which is caused by the Cabibbo-suppressed part of the four-fcrmion operators not included 
in the operator basis ([|). However, the impact of this process in _B — > {K,K*) transitions is negligible Q. 
The amplitudes of meson decays are induced by the quark transition 6 — s- s through the vector = 37^6, 
axial-vector = sj^^b, tensor T^,^ = sa^,yb, and pseudotensor = qa^^^^h currents, with the following 
covariant structure p5| 

< P{M2,P2)\V^mP{Mi,Pi) > = f+{q^)P^ + /-(g')g^, 

< V{M2,P2,e)\A^mP{Mi,Pi) > = ie*''[f{q^)gt.a + a+{q^)p,„P^ + a^{q^)pi^q^], 

< P(M2,P2)|T^.(0)|P(Mi,pi) > = ~2i s{q^) (pi^P2. -Pi.P2m), 

< y(Af2,P2,e)|r^.(0)|P(Mi,pi) > = ie*" [g+{q^)e^,,^f3P^ + g-{q^)e^,,^pq^ + go(q^)p^^e^,,f3^p'lpl], 

< P{M2,P2)\TlM\P{MuPi) > = s(g2)e^^„^P"/, 

< y(M2,P2,e)|T^',(0)|P(Mi,pi) > = g+{q^){etP^ - e;P,) + g^iq^){e:q^^ e;q,) 

+5o(9^)bie*)(Pii/P2A< (6) 

where q — pi — p2, P ~ pi + p2- We use the following notations: 7^ = ?7"7^7^7'^, a^^ — ^[7^, 7i/], e^^^'^ = — 1, 

75^^. = -5ep.a/3fT"^ and Spi-f'^rri") = ^te''''"''- 

Another frequently used set of the form factors is connected with the set (||) as follows 

F^{q') = U{q'), 

Foiq^) = U{q') + qy-{q^)/{Pq), 
FT{q^)^^{M,+M2)siq^), 
Viq^) = (Ml + M2)g{q^), 
A,iq') = f{q^)/{A'h+M2), 
A2{q^)^~{Ah+M2)a+{q^), 

Aoiq') = [q'a-{q') + f{q') + {Pq)a+{q^)]/2M2. 
Tl{q') = -g+{q^)l2, 

T2{q') = -{g+{q') + q'g-iq')/iPq))/'2, 

Tsiq') - (Ml + M2)2[g_(g2)/(Pg) _ h{q')/2]/2 (7) 

The relativistic invariant form factors encode the dynamical information about the decay process and should 
be considered within a nonperturbative approach. We investigate the meson form factors within a dispersion 
formulation of the relativistic constituent quark model (QM) model [ pT[ |. 

Let us consider the transition from the initial meson 5(^2)9(7713) with mass Mi to the final meson (7(7711)5(7713) 
with mass M2, induced by the quark transition 7772 —>■ nii through the current q{mi) J ^i^^-^q[m2) ■ For the 
transition Bu — > {K,K*) one has 7712 = mt, mi = rus and 7773 = niu- The CQ structure of the initial and 
final mesons is described by the vertices Fi and F2, respectively. The initial _B-meson vertex has the spinorial 
structure Fi = 775 Gi/y/Nc, where Nc is the number of colours; the final meson vertex has the structure 
F2 = ij5 G2/\//Vc for a pseudoscalar state and F2;^ ~ [A-f^ + B{ki — fcs)^] G'2/-\//Vc, with A = —1 and 
B = 1/ (\/s2 + ^771 -f 7773) for an 5*- wave vector meson. At < the spectral representations of the form factors 
have the form |^ 



16^2 J J --^^^-^AV2(si,s2,g2)' 

where the wave function Lpi{si) ~ Gi{si) / {si — Mf) and 

_|_ 2\ S2{m\ + 7772 ^ q^) + q^{fn\ + 7773) — (77I1 - 77l2)(777,i — 77l|) A^/^(S2, W^, 777i ) A^/^ ((J^ , ?T7i, TOj) 

^'^''^ 2^1 2^1 
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and A(si, S2, S3) = (si + S2 — S3)^ — 4siS2 is the triangle function. 

The Eq. accounts for only two-particle singularities in the Feynman graphs. For ground-state pseudoscalar 
and vector mesons built up of constituent quarks with the masses niq and m^, the function (p{s) can be written 
as 

, TT \/^"-K-"^|)"u;(fc^) 
where k = \^^'^{s,m^,m'^)/2y/s and w{k^) is the ground-state S-wave radial wave function, normalized as 

dke\wie)\^ = 1. 

The unsubtracted double spectral densities /i(si, S2, g^) of the form factors read [^ : 

s = 2 [mia2 + m2ax + 7773(1 — ai — Q!2)], (10) 

/+ + /_ = /l = 27771S -I- 4a2[s2 - (7771 - 7773)^] - 2r773S, (11) 

/+-/_ = /2 = 27772S + 4ai[si - (7772 - 7773)^] - 27773s, (12) 

g = -As - ABP, (13) 

~9+ + ~g- = Ah - 8/3 + 8B(777i + 7^3)/?, (14) 

9+ - 9- = Ah + 85(7772 - 7773)/?, (15) 

a+_D - a_£i = -2.S + 4i?C2ai + Q!i2Co, (16) 

a+_D + = -4A (27773 + BCi)ai + anCo, (17) 

fa = -4yl[?77 1 7772 7773 -|- ^(^2 - mj - ml) + ^(si - 7772 - ^l) ^ - m\ - ml)] + Co/3, (18) 

gOD = -8Aai2 - SB [-7773Q!l + (7773 - 7772)011 + (7773 + 777l)ai2], (19) 

ai = [{si + S2- S3){s2 - m\ + 7773) - 2s2(si - 7773 -I- 7773)] /A(si,S2, S3), (20) 

a2 = [(si + S2- sz)[si - 7772 + "^s) ^ 2si(s2 - m\ + ml)] /A(si,S2, S3), (21) 

/3 = i [27773 ^ Ctl{si - 7772 + "^3) ^ "2(52 - 777^ + 7773)] , (22) 

an ^ a\+ 4/3s2/A(si, S2, S3), ai2 = aia2 - 2/3(si + S2 - S3)/A(si, S2, S3), (23) 

Co = -8yl (7772 - 7773 )+4BC3, Cl = S2 - (7771 -I- 7773)^ (24) 

C2 = Si - (7772 - 7773)^, C3 = S3 - (7771 -I- 7772)^ - Ci - C2. 

We label with a subscript 'D' the double spectral densities of the form factors which require subtractions. 
The subtraction procedure has been fixed by matching the Xjmq expansion of the form factors in the quark 

model to the corresponding expansion in QCD in leading and next-to-leading orders for the case of a meson 
transition caused by heavy-to-heavy quark transition ||22| . 

The double spectral densities with properly defined subtraction terms read 

/ = /d + ml - si) + {Ml - S2)]g, (25) 

. , Vir+ ( Ml-s^ Mi-S2 \ g 

a+ = a+D + , _ , — ^72 1=— + 1=— t;^ (26) 

. yi^ - yil f Ml - si , Mj -S2 \9 

a- = a^D + , ^ — + o' (27) 



(V^+x/iJ)2-S3 V x/^ y 2' 

1 fM?-si M|-S2\. 

30 = 90D + , ^ , V + V ff- (28) 



where 



The structure of the HQ expansion in LO and NLO 1 / 77iQ orders of the form factors given by the Eq. (||) 
with the spectral densities (p^refgO) agrees with the corresponding structure of the HQ expansion in QCD 
p4|| , provided that the functions <~p{si) are localized near the qq threshold with a width of the order Kqcd p2t . 
Moreover, for the case of meson decays induced by a heavy-to-light quark transition the dispersion formulation 
provides the form facotrs which satisfy the leading-order relations between the form factors of the vector and 
tensor currents near zero recoil given in [25|] . 

As the analytical continuation to the time-like region > is performed, in addition to the normal contri- 
bution which is just the expression (H) taken at q^ > 0, the anomalous contribution, described explicitly in [|l] 
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emerges. The normal contribution dominates the form factors at smaU and vanishes when = (7712 — ^i)^, 
while the anomalous contribution is negligible at small > and steeply rises as q^ — > (rn2 — mi)^. 

Notice that since the dispersion quark model is based on taking into account only two-particle qq intermediate 
states in the amplitude of the interaction of the qq constituent quark pair with the external field it is conceptually 
close to the light-cone quark model LCQM |8j. In particular, the form factors of the LCQM |^ can be rewritten 
at < as double spectral representations similar to the dispersion model. One finds that at < the form 
factors which are given by the unsubtracted spectral representations in the dispersion formulation are the same 
as in the LCQM. At the same time, the LCQM form factors /, oi, 02, and h are different from the dispersion 
quark model form factors and do not develop a correct heavy quark expansion in the next-to-leading I/uiq 
order. 

For evaluating the form factors we need to specify the quark model parameters such as the constituent quark 
masses and the wave functions. In Ref. | [Tl[ | we have run calculations of the mesonic form factors (^) adopting 
different QMs for the radial wave function u'(fc^) appearing in Eq. (|^), in particular: a simple Gaussian ansatz 
of the ISGW2 model ||2^ and the variational solution of the effective qq semi-relativistic Hamiltonian of 
Godfrey and Isgur (GI) |28(|. These two models differ both in the shape of the radial wave function, particularly 
at high momenta fc, and in the values of the quark masses (see Ref. The results of our calculations 

have shown that the mesonic form factors (|^) are sensitive both to the high-momentum tail of the meson wave 
function and to the values adopted for the quark masses (see also Refs. p9| , |30| ). 

In order to obtain more reliable predictions for the form factors, we require the QM parameters to be adjusted 
in such a way that the calculated form factors at large are compatible with the lattice results [ p^p^ . We 
have found that the best agreement with the lattice data at large q^ is obtained for the quark masses and wave 
functions of the GI model with a switchcd-off one-gluon exchange (GI-OGE). The constituent quark masses and 
the average momenta squared characterizing the GI-OGE model are given in Table |. Table ^ presents simple 
fits to the calculated form factors. 

Fig. |l| shows the GI-OGE form factors versus the available lattice data for B K* p^ , p^ . For comparison, 
we also present the form factors obtained within the ISGW2 exponential ansatz for the soft meson wave functions 

0- 

A good agreement of our QM predictions obtained with the GI-OGE wave functions with the results of lattice 
simulations is not surprising: the strong long-distance physics is dominated by the confinement mechanism, and 
therefore it seems quite natural, that soft wave functions which take into account the effects of the confinement 
scale provide the form factors in agreement with the lattice QCD results. 



TABLE I. Constituent quark masses (in GeV) and the average momentum squared (in GeV^). 



Ref. 


niu 


ms 


rric 


mt 


< >K 


< P >K* 


< >B^ 


GLOGE (27) 


0.22 


0.42 


1.65 


5.0 


0.17 


0.17 


0.26 



TABLE II. Parameters of the fit fi{q'^) = /i(0)/[l — uig^ -I- a2q'^] to the B {K,K*) transition form factors in the 
GLOGE model. 



Decay 


B ^ K 


B ~* K* 1 




/+(0) 


/-(O) 


s{Q) 


3(0) 


/(O) 


a+(0) 


a_(0) 


HO) 


9+iO) 


ff-(O) 


Ref. 










(71 


(71 


(71 


(71 


(71 


(71 




(72 


(72 


(72 


(72 


(72 


(72 


(72 


(72 


(72 


(72 


GLOGE 


0.33 


-0.27 


0.057 


0.063 


2.01 


-0.0454 


0.053 


0.0056 


-0.3540 


0.313 




0.0519 


0.0524 


0.0517 


0.0523 


0.0212 


0.039 


0.044 


0.0657 


0.0523 


0.053 




0.00065 


0.00066 


0.00064 


0.00066 


0.00009 


0.00004 


0.00023 


0.0010 


0.0007 


0.00067 



7 




Recently, a lattice-constrained_parametrization for the B — > K* form factors has been developed in the whole 
range of accessible values of q-^ . It is based on the Stech parametrization of the form factors obtained within 
the constituent quark picture, HQS scaling relations near — q^^^x ^iid a single-pole behavior of Ai suggested 
by mg-scaling relations at = from the LCSRs in the HQ limit. Parameters of the single-pole fit to the form 
factor v4i(0) = 0.291q q3 and Mi = 6.8^q 4 GeV are found from the least-x^ fit to the lattice QCD simulation in 
a limited region at high values of q^ . Such a parametrization, though still phenomenological, is also consistent 
with the dispersive bounds of Ref. jl^ and therefore obeys all known theoretical constraints. It should be taken 
into account however that the lattice-constrained parametrization is an approximation: in particular, it suggests 
the relation Ti{q^) = (1 — q^ / Pq)T2{q'^) which can be also translated into g+{q'^) = —g-{q^)- In dynamical 
calculation within QMs or LCSRs these relations are fulfilled within 10% accuracy but are never found to be 
exact. Nevertheless, approximate Stech's relations combined with a monopole fit to Ai exhibit surprisingly 
good agreement with the lattice points at large q^ (see Fig. |]). 

The form factor Ai calculated in our approach for the GI-OGE model wave functions is found to have a 
behavior very close to the single-pole function with the parameters ^i(O) = 0.326 and Mi = 6.86 GeV in 
agreement with an assumption of Ref. p^ ] . The results on most of the form factors are within 5% agreement 
with the parametrizations [0 except for the form factor V which turns out to be at zero recoil some 15% 
smaller in our calculations. 

Table [II compares predictions on the form factors from various approaches. One can see that our results 
agree with those of the LCSR of Ref. jlj] . The form factors of another version of the LCSR Ml have different 
behavior which disagrees also with the lattice results at large q^: namely, at q^ ~ 16 GeT^^the form factors 
Ti, T2, and Ai turn out to be considerably larger than the lattice points, and the form factor ^0 to be too small. 
The form factors T2 and Ai in the 3ptSR approach |jl^ are decreasing with q^ in contradiction with the results 
of other approaches and lattice simulations. 

Let us notice, that an approximate relation g+ = — g_ found to describe well the lattice points at large q^ 
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and extended in to the whole kincmatically accessible region, might signal that relations motivated by the 
heavy-quark symmetry also work with a reasonable accuracy in the B ^ K, K* case. In fact, the HQ expansion 
of the form factors g+ and 5- gives 5+ = ^^±^£iw{oj)[l + 0{l/mQ] and .g_ = -Mi^^j^{u:)[l + 0{l/mQ\. 



2VM1M2 

An approximate relation (7+ — —g^ can be obtained from this expansion in the limit <^ M b only if the 
generically different combinations + 0(1/toq] in and evolve in this limit to the same function 

S^B^K* which however goes far from the Isgur-Wise function. Let us assume the leading-order IW relations for 
the form factors with the IW function replaced by the function £,b^k and £,b^K' for £? ^ X and B — > K* 
transitions, respectively. The process-dependent functions £,b^k,k* determined from the GI-OGE QM results 
for Ti and Fi through the relations 



Mb + M 



K 



Mf 



M 



(29) 



K 



are shown in Fig. g. The deviations for other form factors of i? ^ K* transition found through the LO HQS 
relations with S^b^K' from the lattice-constrained parametrizations can be as much as 20%. 

Summing up, the dispersion quark model calculates the form factors in the whole kincmatically accessible 
decay region. The form factors of the dispersion quark model have the following properties: they develop a 
correct HQ expansion in the leading and next-to-leading 1/toq orders in agreement with QCD in heavy-to- 
heavy transitions provided the soft wave function is concentrated in the region of the confinement scale; for the 
case of heavy-to- light transition they have the correct scaling properties at small recoil and obey the LO 1 /mg 
relations between the form factors of V, A, and T currents pq] ; and numerically they agree with the lattice 
results at large for the B K* transitions. Thus we expect the dispersion quark model form factors to be 
reliable in the whole kincmatically accessible region. 

In the next sections we use the QM form factors evaluated with the GI-OGE wave functions and the lattice- 
constrained parametrizations of Ref. for analyzing the decay rates and asymmetries in rare B decays. 



TABLE III. Comparison of the results of different approaches on the form factors Ti, T2, Ai and Ao at 
f = 16.2 GeV^. 



and 



Ref. 


Ti 


T2 


Ai 


Ao 






5^ = 


= 16.2 GeV 


q' = 16.2 GeV 


g^ = 


= 16.2 GeV'^ 


q'=0 


q' = 16.2 GeV" 


LCQM 






0.155 


0.53 


0.26 


0.26 


0.45 


0.32 




3ptSR 


12 


1 


0.19 ±0.03 


0.5 ±0.05 


0.13 ±0.03 


0.37 ±0.03 


0.23 ±0.03 


0.3 ±0.03 


1.0 ±0.05 


LCSR 


13 




0.18 


0.84 


0.35 


0.36 


0.65 


0.27 


0.64 


LCSR+Ls 


it 


14 


0.15 ±0.04 


0.54 


0.22 










Lat+ 








0.57 ±0.04 


0.21 ±0.04 


0.29 


0.44 ± 0.04 


0.33 


1.28 ±0.07 


GI-OGE 


0.177 


0.53 


0.248 


0.33 


0.44 


0.44 


1.20 


ISGW2 


0.142 


0.46 


0.23 


0.26 


0.42 


0.35 


1.08 
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IV. DIFFERENTIAL DECAY RATES AND LEPTON ASYMMETRIES 



In this Section we present formulas for the differential decay rates, forward-backward asymmetries and lepton 
polarization asymmetries obtained for both nis ^ and mg ^ for the transition induced by the effective 
Hamiltonian Eq. (|) in the case B {K,K*){e+£-) and by Eq. (|) in the case B {K,K*){vv). 

Our formulas for the differential decay rates and forward-backward asymmetry coincide with the corresponding 
formulas of Ref. and reproduce the formulas of Refs. [|lO[^ in the case = and m; ^ and those of 
Ref. p^ ] in the case rus — 0,mi — 0. For lepton polarization asymmetries our expressions in the limit — 
coincide with the results of [lOlJl^ . 

Introducing the dimensionless kinematical variables s = q^/M^ and t = {Pb — pi+)'^ /Mg, the double differ- 
ential decay width for the rare decay B K£~^£~ can be cast into the form 



where 



d^r{B K£+n _ G%Ml\V,lVtb\'al^ 



dsdt 



2567r5 



-n/3p + 2rhSp 



(30) 



/3p = 



C;(/{mi„q')f+{q') + 2{mb + ms)CT^{mi,)s{q') + \CioAimb)f+{q')\ 



11= (i- 1) (t - f) + It + m(l + f + m-s-2i), 



SP = \C,oa\' 1 + ^ - o + (1 - r)Re[f+{q')fliq')] + i;\f-i<l 



(31) 



with r = {Mk/MbY' and m = {ttii/Mb)'^ ■ After integrating over t from tmin = [1 + r -f 2rh — s — 
v'l-4m/sAi/2(l, s, r)]/2 to i^ax = [l + f + 2rn-s + ^\ - Am/ sX^''^{\, s, f)]/2, where A(l, s,f) = l + r^ + 



S2 



2f — 2s — 2fs, one obtains the invariant dilepton mass distribution 



ds 15367r5 ^ ' V , , ; 

In the case of the decay B K*£^£~ one has 



2m 



1 + — ) X{l,s,f)(3p + I2mdp 



d^r{B ^ K*e+£-) Gj,Mlj\VtlVtb\^a^ 



dsdi 



5127r5 



with 



is + 2m)X{l,s,f)+2m |G(g^)r + 



S + 2771 - 



n 

2f 



2\\2 



\Fiq')\ 



XHl,s,f) 



J — 1 + f , 



n|H+(g2)|2 + : ^URiq^l 



2f r 

,2 



\Giq')\' 
\F{q')\' 

R{q^) = Re 



2s[2t + s - f - 1 - 2m]Ri{q^) 

Cg{/{mb,q jMpgiq ) '- — g+{q^) 



+ \CioA{fnb)MB9{q^)\' 



r«^fff 2J{'f') 2Cj^{mb) mb ~ ms --.u f 2-, 



CwAinib) 



Clif{mb^e)MBa,ie) - ?^!2^!I^-^ . ..2. 



C^ffim. o^^Zili _ 2C7^(m,)mfc-m o .^2^ 



9V 

2Cj^{mb) mb — m 
i Mb 



Mb 

+ \CioAimb)MBa+iq^)\^, 

Cl(/{mb,q^)MBa+{q^) 



■B+{q' 



+ \CioA{fnb)\'Re[a+{q')r{q% 



Ri{q^) = Re 



r«^ffi '2\AT ^ 2n 2C7^{mb)mb + ms 2^ 

Cgi/(mb,q )MB9{q ) '- — g+{q ] 

s M R 



CioAimb) 



-Re 



Cgi/{mb,q )— '- — (l-r)Bo(q ) 

i\Ib s Mb 



Mb 

[CwA{mb)MBg{q^)\* 



(32) 



(33) 
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B+{q^) = -sM: 
Sv 



r2 



i-r 

2\ 



\c- 



X{l,s,f){-2\g{q')MBV 



A(l,s,f) 



2{l + k)-s 
Af 



^|a_(<z2)Af5|2 + _Li?e[/(g2)a;(g2) + J{q^)a*_{q^)] + l^Re[MBa+{q^)MBa*_{q^)] \ . (34) 



where now f = {Mk* /MbY' ■ After integrating over the variable t we find 



ds 



15367r5 



1 + 



2m 



(35) 



where 



/3y = 2A(l,s,f)J|G(g2)|2 
A(l,s,f 



2s 



(1 - f - s)2 
4f 



|i^(g^ 



2\|2 



A^(l,^,r) 
4f 



H+{q 



2\\2 



2f 



is-l + f)R{q^). 



(36) 



The effective Hamiltonian for the B (X, iir*)(i/z/) transition (||) may be obtained from the corresponding 
Hamiltonian for the B — > {K,K*){i~^£~) transition (Q) by the following replacements: 



m ^ 0, Cr-y 0, q 



9V 



X(xt) 
sin^(0w)' 



^(a:t) 
sin^ {9w) 



(37) 



Hence, expressions for the decay rates in i? ^ (K, K*){i'iy) can be obtained from the corresponding formulas 
for B ^ {K,K*){i+i~) by the replacement (|^). 

For the decays B K*£^£~ a very interesting quantity is the forward-backward (FB) charge asymmetry 
^_fb(s), which is defined as 



ApBis) 



1 



dT{B K*i+£-)/ds 



\ d^r{B ^ K*£+£-) ? ^^ d^r(B^K*£+£- 
dcos(O) — — — / acos(y) 



dsdcos{6) 



dsdco8{9) 



(38) 



where 9 is the angle between the charged lepton £^ and the _B-meson directions in the rest frame of the lepton 
pair. As is well known, the FB asymmetry is sensitive to the parity structure of the electroweak interaction. At 
low values of q^ the parity-conserving photon exchange is expected to dominate and therefore the FB asymmetry 
should be small; on the contary, at large q^ the contribution of the parity- violating Z- and VF-boson exchanges 
becomes relevant, leading to a large asymmetry. Moreover, the FB asymmetry is sensitive to the relative sign 
of the Wilson coefficients ||] and therefore its measurement could be used as a probe of the new physics beyond 
the Standard Model. Explicitly, one has 



Afb{s) 



3s^l -4m/sAi/2(i^s^f)i?i(q2) 



(1 + 2f ) /3y + 12m,5y 

Finally, we will consider also the longitudinal lepton polarization asymmetry Pl(s) defined as 



dV/ds 



dT{h^-l) dT{h^+l) 



ds 



ds 



(39) 



(40) 



where h = +1(— 1) means right (left) handed charged lepton £ in the final state. In case of the rare decay 
B K£+£- one has 



2y/l-4rh/sX{l,s,f) 

1 + 2m\ ^-^^ ^ 127h5p 



Re 



{ [c;(/{mi„ g')/+(<?') + 2(mf, + m,)CTy{mt,)s{q^) 



(41) 



whereas for the process B — > K*£^£ one gets 
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2y/l ~ Am/s 



where 



Raiq^) = Re 
Rpiq^) = Re 
RH^{q^) = Re 
Riiiq^) - Re 
Re 



, A2(l,s,f) 



2A(1, s, f)sRG{q^) + ( 2s + \. ) RF{q^) 



4f 



4f 



A(l,g,r) 
4f 



{s-l + f)RR{q') 



Cgi/{mb,q )MB9{q ) — " 



M 



B 



-9+iq') 



C,i.'{mb,q)—- ^ ) 



IB 



Mf 



-B+iq') 



Cgl/(mb,q )MBa+{q )- ' d ^„ ^ 

Cl{/{mb,q^y 



f{q^) 2C'7y{mb) TUb - rus 



Mb 



Mb 



il-r)Bo{q') 



Cgyimb.q )MBa+[q ) '- — ° ^ 



Mf 



'-B+iq') 



[CioA{mb)MBg{q^)] 

CioA[mb) 
[CioA(™fc)M_Ba+(g^)]* 

[CioA(mfc)MBa+(g^)] 
fiq'-'* 



(42) 



Cioyi(m6)- 



Ms 



(43) 



Notice that for a left-handed massless neutrino Eq. (40) yields = —1. The same result can be obtained also 
fromEqs. (p7|M P). 



V. NUMERICAL ANALYSIS 

In this section we analyze the decay rates and lepton asymmetries in the Standard Model: the Wilson 
coefficients at the scale /i ~ rUb are given in Section 2, and the B K,K* transition form factors of the 
dispersion quark model evaluated with the GI-OGE model wave functions are adopted. For the B K* 
transitions we also consider the lattice-constrained parametrizations of the form factors of Ref. . We denote 
the two sets of predictions as QM and Lat, respectively. 



A. Decay rates and dilepton distributions 



First, let us evaluate the \Vts\ 
§ we find 



Combining our QM prediction for the r2(0) with the CLEO data on S K*j 



0.038 ± 0.005'=''P. 



(44) 



A similar analysis with the lattice-constrained parametrization for T2 jlj] yields \ Vts\ = 0.041±0.005*'^±0.005™p. 

The predictions for the dilepton distribution in B ^ {K, K*)£~^£~ decays are reported in Fig. ^, where the 
non-resonant contributions are also shown separately. The total decay rates turn out to be at least one order of 
magnitude larger than the non-resonant decay rates. However, the resonant contributions are strongly peaked 
in narrow regions around their masses, so that outside these regions the resonance influence is almost negligible. 
This fact allows one to reliably separate the resonant contribution from the non-resonant one, which is mostly 
interesting as it contains the information on the Wilson coefficients and in principle allows one to measure these 
coefficients. The resonance phases are chosen in accordance with the analysis of Ref. 

Table IV summarizes our predictions for the non- resonant branching ratios. The branching ratios obtained 



with the QM and Lat sets of the form factors are given in units |Vts/0.038p and |Vts/0.041p, respectively, such 
that the B — > K*j branching ratio evaluated with Ct^{ii = 5 GeV) in each case is normalized to the central 
CLEO value 0] if this factor is equal to unity. 

Note that the transitions B K*fj,^fi~ and B K*e^e^ have different rates, because the amplitude 
B K*l^l'^ has a kinematical pole at = 0, which makes the corresponding decay rate very sensitive to the 
lower boundary of the phase space volume (g^ = 4m^), while the amplitude B K£'^£~ is regular at = 
and, therefore, insensitive to the mass of the light lepton. 

One observes a strong sensitivity of the differential decay rates in B ^ K* transitions at low q^: the decay 
rates are results of the interference of various form factors and thus are sensitive to the details of their q^- 
behavior. One can see that the Lat and QM form factor sets which both provide a reasonable agreement with 
the lattice results at large 5^ , yield sizeable deviations in the differential distributions at low q^ and hence in the 
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branching ratios. Thus for deriving accurate predictions for the branching ratios one needs accurate knowledge 
of the form factors at fow . 

The decay rates of the B ^ K modes turn out to be more stable with respect to variations of the relevant 
form factors (cf. [pl| ) and thus might be more perspective for extracting Vts from rare semileptonic decays. 




0.2 OA 0.6 O.S 0.2 0.4 0.6 0.8 




FIG. 3. Differential decay rates lO^d Br/d{q^/Ml) in the decays B -> K*l+r: &. B K/i+fj.". h. B ^ K* fi+ fi' . 
c. B ^ Kt+t- , d. B ~* K*r+T'. Solid - QM form factors (GI-OGE), for \Vts\ = 0.038; dotted - lattice-constrained 
parametrization of for \Vts \ = 0.041. Thick lines - nonresonant parts, thin lines - total. 



TABLE IV. Non-resonant branching fractions of the radiative and rare semileptonic B decays. The branching 



fractions are evaluated using Tie// at the scale fi — 5 GeV. Uncertainties connected with a relevant choice of this 
low-energy scale fi ~ mj, are not shown. 



Decay 
mode 


QM 

x\Vts/omsf 


Lat 
x|Vt,/0.0411^ 


i 

x|Vf,/0.033p 


Exp. 


B K*j 


4.2 X 10"'' 


4.2 x 10"' 


(4.9 ± 2.0) X 10"" 


(4.2 ± 1.0) X 10" 


- b 


|l| 


B K*e+e- 


1.50 X 10"" 


1.45 X 10"" 


(2.3 ±0.9) X 10"" 


< 1.6 X 10"" 1 


52 




B K'li+fi- 


1.15 X 10"" 


1.1 X 10~" 


(1.5 ±0.6) X 10"" 


< 2.5 X 10"" i 






B ^ K't+t- 


1.0 X 10~' 


1.1 X 10"' 






B ^ K*J2 ^i^i 


1.5 X 10"-' 


1.4 X 10"" 


(1.1 ±0.55) X 10"" 








B Ki+r 


4.4 X 10" ' 




(4.0 ± 1.5) X 10" ' 


< 0.9 X 10"" [32 


B ^ Kt+t- 


1.0 X 10"' 








B ^ UiUi 


5.6 X 10"" 




(3.2 ± 1.6) X 10"" 
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B. Forward-backward asymmetry 



The forward-backward dilepton asymmetries in B ^ K*i'^£^ decays are presented in Fig. For comparison, 
we show also nonresonant Afb evaluated with the form factors for the ISGW2 quark model parameters, and 
obtained assuming the HQS relations between the form factors. One can observe a strong sensitivity of the 
asymmetries to the specific details of the behavior of the relevant form factors. Notice that the maximum of 
the asymmetry at q^/Mg ~ 0.1 is mainly proportional to the ratio (V/Ai)^. That is why the Lat Afb turns 
out to be larger in the maximum than the QM asymmetry. Nevertheless the general trend of the behavior of 
Afb for all considered sets of the form factors is similar: the nonresonant asymmetry is positive at low g^, has 
a zero at g^/Af^ ~ 0.15, and then becomes negative irrespective to the details of the form factor behavior. Let 
us point out that maximum absolute value of the (negative) asymmetry is attained at g^/Af^ ~ 0.62 where 
|Afb| — 0.4 both for the QM and Lat sets of the form factors. This is considerably smaller than |^_f_b| ~ 0.6 
reported in [Q. This difference is traced back to a very specific behavior of the form factor T3 in |]l^ which 
contradicts to the results of other approaches and to the approximate HQS relations between the form factors 
(see also discussion in [0). 




FIG. 4. Forward-backward asymmetries in B ^ K*t^i~ transitions, a. B K* fi'^ fi^ {e'^ e~ ) , non-resonant, b. 
The same, total, c. B —* K*t^t^ , non-resonant, d. The same, total. Solid - GI-OGE, dashed - ISGW2 models, dotted 
- lattice-constrained parametrization of |17[|, dash-dotted -HQS relations. 
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C. Lepton polarization asymmetry 

Fig. 1^ shows lepton polarization asymmetries Pl for massless and massive leptons. For understanding the 
behavior of Pl it is important to take into account the relationship between the Wilson coefhcients in the SM: 

C'r^irrib) < CioA{mb) ~ -Cgv(m{,). (45) 

In the case of the transition B K£^£^ , £ — a. simple analysis of Eq. ( |4]] ) yields the following behavior 
of the nonresonant P^: Pl is equal to zero at — Amj and = {Mb — MkY '^^'^ fo kinematical reasons, 
and in the intermediate region of q^, Pl steeply goes down to the value Pl — 2CqvCioa/ {Cly + C^qa) — ~^ 
independently of the particular behavior of the B —> K transition form factors. A weak g^-dependence of the 
nonresonance Pl is due to the function h{mc/mb,q'^ /mf) in Cg(/ . 

In the reaction B — > K*£'^£^, £ — the situation is a bit different: Now the term in Ti-eff proportional 
to a small Cj^ contains a photon pole at = and thus a parity-conserving photon exchange dominates the 
decay at low providing a small value of Pl- At large one finds Pl — —I because of just the same reason 
as in the B ^ K case with the only difference is that the kinematical zero at = {Mb — Mk'Y is absent. In 
the intermidiate region of g^, the nonresonant Pl is an interplay of the parity-conserving and parity- violating 
terms yielding a negative Pl smoothly falling from to —1 in a way largely independent of the particular 
g^-dependence of the transition form factors. 

In the total Pl in the reactions B {K, K*)£'^£^ , £ ^ ii,e the t/j and ip' resonances appear as sharp peaks 
on a smooth nonresonance background. 



B K\l*iL 



o.s 





FIG. 5. Longitudinal lepton polarization asymmetry {Pl) in the decays B — > {K, K*)£'^£~ : a. B — > K/j,'^ {e'^ e~). 
b. -B — > K* fj,^ {e^ e~). c. B — » Kt^t~ . d. B ^ K*t^t~ . Solid - GI-OGE, dotted - lattice-constrained parametriza- 
tion of |17|- Thick lines - nonresonant, thin lines - total. 
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The results of our calculation shown if Fig. ^ as well as the results of Ref . 11 1 correspond to the picture 
described above, whereas Pl reported in |13| has a different behavior with Pj, — —0.6 at large which seems 
to be very doubtful. 

The lepton polarization asymmetry in the case B (K, K*)t^t^ in general follows the trend of the 
light leptons case with an important difference: the nonresonant Pl does not go down to the value ~ — 1 in the 
kinematically accessible region. The P^ again turns out to be largely insensitive to the meson transition form 
factors. In the total Pl one observes only the tp' peak in the kinematically accessible region. 



VI. CONCLUSIONS 



We have analyzed rare semileptonic transitions B {K, K*) within the Standard Model adopting two models 
for the relevant form factors: a relativistic constituent quark model, formulated in a dispersion form, and the 
lattice-constrained parametrization of Ref. [pT| . Our main results are as follows: 

• We have presented a dispersion quark model calculation of the B ^ K, K* transition form factors in 
the whole kinematical range of q^. Adopting the quark masses and the wave functions of the Godfrey- 
Isgur model for the hadron spectrum with a switched-off one-gluon exchange potential for taking into 
account only the impact of the confinement scale, we have found the resulting form factors to be in good 
agreement with the lattice simulations at large q^. 

The form factors in the dispersion quark model develop the correct expansion in the leading and next- 
to-leading l/mg orders for the heavy-to-heavy decays, and satisfy the relations between the form factors 
of the vector, axial-vector, and tensor currents valid in the region near the zero-recoil point in case of 
heavy-to-light decays. In addition, the form factors are compatible with known analytical constraints. 

Hence, the form factors of the dispersion quark model obey all existing rigorous theoretical constraints 
and agree nicely with the results of lattice simulations for the B — > K* decay at large q^. Moreover, the 
dispersion quark model form factors for the B — > K* transition agree favorably in the whole range of 
< q^ < {Mb — Af^)^ with a lattice-constrained fit [0 based on the constituent quark picture |^ and 
an assumption on a single-pole behavior of ^i(g^). Thus we expect to have reliable form factors in the 
whole kinematically accessible decay region. 

• We have performed a detailed analysis of the non-resonant decay rates and asymmetries in _B — > 
{K,K*) {t^i~,vv) decays in the Standard Model and obtained predictions for all exclusive channels 
using our GI-OGE form factors and the lattice-constrained fit to the form factors oi B ^ K* transition. 

i. Combining our QM result for T2(0) with the central CLEO value for B K*^ [Q we estimate the 
central value \ Vts\ — 0.038. With the lattice-constrained parametrization of the form factors ||l^] one finds 
the central value \Vts\ = 0.041. 

ii. The results of the non-resonant branching fractions, obtained within the two sets of the form 
factors, are in good agreement if the relevant |Vts| is used in each case. Nevertheless, a better knowledge 
of the relevant form factors around = is still required. 

iii. The differential dilepton distributions in B K£^£^ decays are less sensitive to the details of 
g^-behavior of the form factors than the corresponding distributions in S — + K*t^£~ processes. Thus, 
the reaction B — > Kfi'^^~ seems to be the most appropriate one for the determination of Vts from rare 
exclusive semileptonic decays. 

iv. The shape of the forward-backward asymmetry in _B ^ K*i.i^^i^ within the SM is almost 
independent of the long-distance contributions: Aps is positive at small q^, has a zero at q^ ~ 0.15M^ 
and then becomes negative at larger q^. On the other hand, the values of Aps in the maximum and the 
minimum are determined by the ratios of the form factors (see also discussion in p^). 

V. The longitudinal lepton polarization asymmetry Pl{B Kfi^ at all kinematically accessible 
g^, except for the end-points and narrow regions near ^ and "0', as well as Pl{B — > K* ji^ at large 
do not depend on the long-distance contributions. In particular, Pl{B K^^^^) ~ 2Cg\/Cioyi/(C|y -I- 
C^Q^)(~ —1 in the SM), and hence Pl directly measures the ratio of the Wilson coefficients Cgy/CioA at 
the sacale /x ~ rnb- Thus, the experimental study of the forward-backward asymmetry and the longitudinal 
lepton polarization asymmetry provides an effective test of the Standard Model and its possible extentions. 

The presented results for the decay rates are essentially based on the lattice-constrained constituent quark 
picture. Further progress in obtaining more accurate predictions by combining these approaches may be ex- 
pected on the following way: with increasing the accuracy of the lattice predictions one can put forward the 
determination of the meson wave functions from the least-x^ fit to the lattice results at small recoils, basing on 
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our proposed spectral representations for the form factors. Then such lattice-constrained quark model would 
provide reliable and accurate form factors at all kinematically accessible . 
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